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ABSTRACT 

We present a new spectrum of the quasar HS 1603+3820 taken 1.28 years (0.36 years in the quasar 
rest frame) after a previous observation with Subaru+HDS. The new spectrum enables us to search for 
time variability as an identifier of intrinsic narrow absorption lines (NALs). This quasar shows a rich 
complex of C iv NALs within 60,000 km s~ 4 of the emission redshift. Based on covering factor analysis, 
Misawa et al. found that the C iv NAL system at z a t, s = 2.42-2.45 (System A, at a shift velocity of 
v s h= 8,300-10,600 km s _1 relative to the quasar) was intrinsic to the quasar. With our new spectrum, 
we perform time variability analysis as well as covering factor analysis to separate intrinsic NALs from 
intervening NALs for 8 C iv systems. Only System A, which was identified as an intrinsic system in the 
earlier paper by Misawa et al., shows a strong variation in line strength (W b s ~ 10.4 A — > 19.1 A). We 
speculate that a broad absorption line (BAL) could be forming in this quasar (i.e., many narrower lines 
will blend together to make a BAL profile). We illustrate the plausibility of this suggestion with the 
help of a simulation in which we vary the column densities and covering factors of the NAL complex. 
Under the assumption that a change of ionization state causes the variability, a lower limit can be placed 
on the electron density (n e >3 x 10 4 cm -3 ) and an upper limit on the distance from the continuum 
source (r < 6 kpc). On the other hand, if the motion of clumpy gas causes the variability (a more 
likely scenario), the crossing velocity and the distance from the continuum source are estimated to be 
Across > 8,000 km s _1 and r < 3 pc. In this case, the absorber does not intercept any flux from the 
broad emission line region, but only flux from the UV continuum source. If we adopt the dynamical 
model of Murray et al., we can obtain a much more strict constraint on the distance of the gas parcel 
from the continuum source, r < 0.2 pc. 

Subject headings: quasars: absorption lines - quasars: individual (HS 1603+3820) 



1. INTRODUCTION 

QSO absorption lines are usually classified by the width 
of their profiles: broad absorption lines (BALs; FWHM 
> 2,000 km s _1 ), narrow absorption lines (NALs; FWHM 
< 200-300 km s _1 ), and mini-BALs (intermediate class 
between BALs and NALs) (Hamann & Sabra 2003). Be- 
cause of their large blueshifts and their wide and smooth 
line profiles, BALs are plausibly associated with accretion- 
disk winds (e.g., Murray & Chiang 1995; Arav et al. 1995; 
Proga, Stone & Kallman 2000). But BALs appear in 
only ~ 10 % of optically selected quasars (Weymann et 
al. 1991). 

On the other hand, NALs can arise from both inter- 
vening absorbers (e.g., cosmologically intervening galax- 
ies, IGM, and quasar host galaxies) and intrinsic ab- 
sorbers (i.e., clouds that are physically associated with 
the quasars). A significant fraction of NALs are thought 
to be intrinsic to the quasars and may arise in gas that 
is outflowing from the central region, based on a statis- 
tical analysis (Richards et al. 1999; Richards 2001). In 
the spectra of radio-quiet and flat radio spectrum quasars, 
intrinsic NALs can appear at a large shift velocities up to 
30,000 km s" 1 (Richards et al. 1999). Because NALs are 
often not saturated, they give a powerful diagnostic of the 
physical conditions of the gas clouds that are physically 



associated with the quasars. 

In order to understand the quasar central region using 
NALs, we must determine which ones are intrinsic to the 
quasars. We can use several criteria to separate intrin- 
sic NALs from intervening NALs; e.g., time variability, 
high electron density, partial coverage, polarization, line 
profile, high ionization state, and high metallicity (Bar- 
low & Sargent 1997; Hamann et al. 1997b, and references 
therein). Among them, two methods, (i) time variability 
analysis and (ii) covering factor analysis, are the most reli- 
able and most frequently applied (e.g., Barlow et al. 1992; 
Hamann, Barlow, & Junkkarinen 1997a; Barlow & Sar- 
gent 1997; Hamann et al. 1997b; Goodrich & Miller 1995; 
Ganguly et al. 1999,2003; Misawa et al. 2003; Narayanan 
et al. 2004; Wise et al. 2004). 

In this paper, we apply the above tests to the absorp- 
tion lines of the bright quasar HS 1603+3820 (z em =2.542, 
B=15.9). First discovered in the Hamburg/CfA Bright 
Quasar Survey (Hagen et al. 1995; Dobrzycki et al. 1996), 
this quasar shows a large number of metal NALs in the 
vicinity of the emission redshift. Dobrzycki, Engels, & Ha- 
gen (1999; hereafter D99) detected 13 C iv NALs at 1.965 
< z a b s < 2.554 in their intermediate resolution spectrum 
(R ~ 3000). The number density of C IV Poisson-systems 4 
with rest-frame equivalent width W res t > 0.15 A is dN/dz 



1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 

2 Department of Astronomy and Astrophysics, Pennsylvania State University, 525 Davey Lab, University Park, PA 16802 

3 Subaru Telescope, National Astronomical Observatory of Japan, Hilo, HI 96720 

4 Groups of C IV NALs that lie within 1,000 km s _1 each other are classified as a Poisson system, because these could be produced in the same 



physical structure (e.g., Sargent, Boksenberg, 



Steidel 1988; Steidel 1990; Misawa et al. 2002). 
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~ 12 at z abs ~ 2.38 in the spectrum of HS 1603+3820 (Mi- 
sawa et al. 2003; M03 hereafter), while the global average 
value is dN/dz = 2.5+og at z abs — 2.40 (Sargent, Boksen- 
berg, & Steidel 1988; Steidel 1990; Misawa et al. 2002). 
Even though the same trends have already been noted sta- 
tistically for the class of radio-loud quasars (Foltz et al. 
1986), such an extreme clustering is significant and prob- 
ably related to the environment around the quasar (e.g., 
intrinsic material, a cluster of galaxies around the quasar, 
or a proximity effect due to the UV flux from the quasar) . 

M03 presented a high-resolution (R/^45,000) spectrum 
of the quasar that covered 9 of the 13 C iv NALs at z a b s > 
2.29, although one of them was apparently not real since 
it was not detected in that higher quality spectrum. Using 
only covering factor analysis, they concluded that the C IV 
NALs at 2.42 < z abs < 2.45 (System A in M03) are intrin- 
sic to the quasar, at velocity shifts 5 from the quasar, v s h = 
8,300-10,600 km s -1 . On the other hand, Systems B, C, 
and D at z abs — 2.48-2.55 could not be demonstrated to 
be intrinsic, although Systems C and D are close (within 
1,000 km s _1 ) to the quasar emission redshift. In fact, one 
of them (System D) is redshifted from the quasar. These 
absorbers (especially the redshifted system) could still be 
intrinsic, but with clouds that happened to cover the emis- 
sion source entirely at the time of observation. Partial 
coverage is a sufficient, but not a necessary, condition to 
demonstrate the intrinsic nature of an absorber. 

In this paper, we present a new spectrum of the quasar 
taken with the same instrument configuration as our pre- 
vious spectrum, 1.28 years after the previous observation 
(0.36 years in the quasar rest frame) with the High Dis- 
persion Spectrograph (HDS; Noguchi et al. 2002) on the 
Subaru Telescope. Our additional spectrum enables us to 
search for time variability of absorption lines as an identi- 
fier of intrinsic NALs, which will not only strengthen our 
classification results, but also constrain several parameters 
of the intrinsic NAL absorbers. 

In §2, we describe the observations and data reduction. 
The methods of line identification and line classification 
are outlined in §3 and §4. In §5, the properties of 8 C iv 
systems, including 4 new systems, are examined in detail, 
and the possible origins of the time variability are dis- 
cussed in §6. We summarize our results in §7. In this 
paper, we use a cosmology with Hq=72 km s _1 Mpc _1 , 
Q TO =0.3, and O A =0.7. 

2. OBSERVATION AND DATA REDUCTION 

We obtained a spectrum of HS 1603+3820 with Sub- 
aru+HDS on 2003 July 7 (UT). We adopt z em = 2.542 as 
the systemic redshift of the quasar, which was estimated 
by matching the spectrum of D99 to the SDSS composite 
spectrum by using the O I A1304, Si n A1307, and C n 
A1335 emission lines (M03). We used a Of! 8 slit width (R 
~ 45,000) and adopted 2 pixel binning along the slit. The 
red grating with a central wavelength of 4900 A was cho- 
sen, which covers the range 3530-4800 A on the blue CCD 
and 4980-6180 A on the red CCD. This configuration cov- 
ers all C IV NALs identified in D99. We took three spectra 
with exposure times of 3,600 s, 3,600 s, and 2,400 s, re- 



spectively. Since the weather conditions were not good at 
the beginning of the observation, we used only the last 
two exposures to produce the final spectrum. We reduced 
the data in a standard manner with the IRAF software 6 . 
Wavelength calibration was carried out using the spectrum 
of a Th-Ar lamp. Because the blaze profile function of each 
echelle order is changing with time, we cannot perform flux 
calibration using the spectrum of a standard star. There- 
fore we directly fitted the continuum, which also includes 
substantial contributions from broad emission lines, with 
a third-order cubic spline function. Around heavily ab- 
sorbed regions, in which direct continuum fitting is diffi- 
cult, we use the interpolation technique introduced in M03. 
We have already confirmed the validity of this technique 
by applying it to a stellar spectrum in M03; the contin- 
uum model is good to 3.3% or better. After binning every 
3 pixels, the signal-to-noise ratio (hereafter, S/N) in our 
6,000 s exposure is <~ 70 per resolution element around 
A = 5450 A. In Fi gure 1, we show the normalized spec- 
trum of the region redward of the Lya emission line of the 
quasar. Echelle order gaps and bad pixels in the detector 
cause the defects seen in the spectrum. The spectrum in 
the region of the Lya forest is not presented because of 
low-reliability of continuum fitting, although some strong 
metal lines were identified by fitting the continuum locally. 

3. LINE DETECTION AND FITTING 

For the present work, we used the line fitting soft- 
ware package minfit (Churchill, 1997; Churchill, Vogt, 
& Charlton, 2003), while we used VPFIT (Webb 1987; 
Carswell et al. 1987) in M03. The minfit code can fit ab- 
sorption profiles using not only redshift (z), column den- 
sity (log TV in cm -2 ), and Doppler parameter (6), but also 
covering factor (C/, described below) as a free parameter. 

D99 detected 13 C iv NALs at 1.88 < z abs < 2.55 in 
their intermediate resolution spectrum. In our new spec- 
trum, we identified 11 of the 13 C iv NALs in that range. 
We also found a new C iv NAL at z abs = 2.2653. We, did 
not detect two C iv NALs found in D99, at z abs = 2.1762 
and 2.5114; the former is unfortunately in the echelle or- 
der gap in our new spectrum, and the latter seems to be 
a false detection because it was noted as uncertain in D99 
and it is not detected in either the spectrum of M03 ( "the 
1st spectrum" , hereafter) or in the spectrum of this paper 
("the 2nd spectrum", hereafter). In total, 12 C iv nALs 
were detected in our new spectrum. 

Following the method of M03, we divided these C iv 
NALs into 8 systems (Systems A-H in Figure 1). Sys- 
tem A contains 5 C iv NALs, while the other systems 
have only one NAL, as identified in D99. These NALs are 
resolved into multiple components in our high resolution 
spectrum (we refer to these as "doublets" in the case of 
C iv, N v, and Si IV lines and "lines" for other single metal 
lines). The first four systems (Systems A-D) had already 
been identified in M03. An additional four systems, at 
Zabs — 1.88-2.26 (Systems E-H), are blueshifted from the 
other systems and were not covered by the 1st spectrum. 
After careful fitting trials with minfit, we found 47 C iv, 



14 Si iv, 6 N v doublets, and 54 other single metal lines 

5 The shift velocity is defined as positive for NALs that are blueshifted from the quasar, which is opposite to the definition given in M03. 

6 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Fig. 1. — The normalized spectrum of HS 1603+3820 (binned by 3 pixels) showing the C IV absorption systems as well as other metal lines. 
The quasar emission lines of Lya, Si IV A1394, and C IV A1548 are marked with downward arrows. The region blueward of Lya emission 
line is not shown, because the spectrum normalization by continuum fitting is not reliable in the Lya forest. The lower line is the la error 
spectrum. The regions affected by either data defects or cchcllc order gaps are masked. 
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associated with the 8 C iv systems. We also fitted four 
systems in the 1st spectrum again, using MINFIT, and con- 
firmed that the results with VPFIT in M03 are consistent 
with the results of MINFIT. 

The fit parameters of these systems are summarized in 
Table 1 for NALs already identified in the 1st spectrum 
and in Table 2 for NALs detected in only the 2nd spec- 
trum. Column (1) is the metal line identification (ID) 
number. Columns (2) and (3) give the observed wave- 
length and redshift. Column (4) gives the velocity shift 
from the quasar systemic redshift, and column (5), the 
Dopplcr parameter. Columns (6) and (7) give the cov- 
ering factors and column densities with la Poisson error 
bars (see discussion in §4.1). Columns (2)-(7) in Table 
1 are the parameters of the components detected in the 
1st spectrum. The corresponding parameters for the 2nd 
spectrum are presented in columns (8)-(13). The same 
parameters for newly identified NALs are listed in Table 
2. Results of a fit for the C iv doublets in System A, car- 
ried out manually (as described in §4.1, below), are given 
in Table 1. Manual fitting was necessary because these 
doublets are so heavily blended with each other and with 
other lines that line fitting with MINFIT cannot be applied. 
For the same reason, the fit parameters for the Si II A1526 
lines in System B are evaluated in the same manner, with 
fitting errors evaluated with the help of the x 2 test. For 
all other doublets, Cf is a free parameter. 

4. IDENTIFICATION OF INTRINSIC NALS 

We describe here two methods of determining that a 
given NAL is intrinsic: (a) covering factor analysis and 
(b) time variability analysis. 

4.1. Covering Factor Analysis 

The covering factor, Cf, is the line-of-sight coverage 
fraction of the absorber over the continuum source and the 
emission line region. The effective covering factor is evalu- 
ated by considering the fraction of photons from the back- 
ground source(s) intercepted by the absorbing gas. The 
normalized flux at wavelength A is given by, 

R(\) = [l-C f (\)]+C f (\)e-TW , (1) 

where r(A) is the optical depth at A. For a resonant dou- 
blet such as C iv AA1548,1551, the ratio of the oscillator 
strengths of blue and red members is 2, as dictated by 
atomic physics. Therefore we can evaluate the covering 
factors as a function of velocity for doublets using the fol- 
lowing equation: 



Cf(v) 



1 + Rl(v) - 2R r {v) 
1 + R b (v) - 2R r (v) 



(2) 



where Rb{v) and R r (v) are the residual fluxes of the blue 
and red members of the doublets in the normalized spec- 
trum (Hamann ct al. 1997b; Barlow & Sargent 1997; Cren- 
shaw et al. 1999). In general, the emission line region and 
continuum source could have different covering factors, as 
discussed by Ganguly et al. (1999). 

To strengthen our results, we evaluate the covering fac- 
tor using the following methods: (i) the pixel-by-pixel 
method - evaluating Cf for each pixel [using equation (2); 
e.g., Ganguly et al. 1999], and (ii) the automatic fitting 
method - fitting absorption profiles using Cf as well as 



log N, b, and z abs as free parameters in minfit (e.g., Gan- 
guly et al. 2003). Unfortunately, neither of these two 
methods is applicable to System A because of heavy self- 
blending (i.e., blending of the blue members of some dou- 
blets with the red members of other doublets). Therefore 
we adopt a third method for this system, manual fitting, 
which we apply only to components 1 and 2, where blend- 
ing is not severe (see Figures 2 and 3). In the manual fit- 
ting method we construct a simple model for the absorber 
with the same free parameters as in the minfit models. 
We then vary the free parameters over a wide range of val- 
ues in small steps, assessing the goodness of the fit first by 
eye, and then by the \ 2 test. 
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Fig. 2.— First: The best-fitting model, with fixed Cf = 1.0, of 
the region around components 1 and 2 of System A in the second- 
epoch spectrum. Model profiles of each C IV doublet and Si II lines 
in System B arc also plotted separately (light solid, dashed, and 
dotted lines respectively). The final model is made to fit the blue 
member of the C IV doublet. The fit of the red doublet member is 
poor, illustrating that a model with full coverage is not acceptable. 
Second: The model that provides the best fit to the second-epoch 
spectrum with Cf treated as a free parameter. This fit is acceptable, 
yielding Cf = 0.45. Third: The model providing the best fit to the 
first-epoch spectrum with the covering factor fixed to Cf = 0.45, 
the value required by the second-epoch spectrum (all other param- 
eters were left free). The model parameters were adjusted so as to 
fit the blue member of the C IV doublet. The fit to the red doublet 
member is poor, however, indicating that the covering factor must 
have varied between the two epochs. Additional details are given in 
§4.1 of the text. Fourth: The model that provides the best fit to 
the first-epoch spectrum with Cf treated as a free parameter. This 
fit is acceptable, yielding Cf = 0.3. 
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Fig. 3. — Velocity plots of the C IV doublets in System A at z a i s 
= 2.4185-2.4455 (v sh = 8, 300-10, 600 km s" 1 ), detected in the 2nd 
spectrum. Vertical dotted lines with ID numbers denote the posi- 
tions of C IV doublets. Smooth lines are models fitted manually, 
using Cj = 0.3 for components 1 and 2, and Cf = 1.0 for the other 
components. The locations of the red and blue members of the 14 
C IV doublets are also marked in the top and the bottom panels to 
illustrate the severe blending. 

The Cf values with their lcr "Poisson" errors, evalu- 
ated by the automatic fitting method, except for the C iv 
doublets in System A, are listed in Tables 1 and 2. The 
Poisson errors are uncertainties in the model parameters 
resulting from only Poisson noise in the spectra. There 
are additional sources of uncertainty, of course, which we 
discuss in detail below. In §5 we present the spectra and 
discuss the properties of the 8 C iv systems to which we 
have applied our analysis. In the case of System A, we 
had to fit the profiles of components 1 and 2 manually. 
We estimated the Poisson uncertainties for components 1 
and 2 by scanning a range of values about the best fit 
for all interesting parameters and looking for a significant 
increase in x 2 (see Lampton, Margon, & Bowyer 1976). 
The uncertainties, which are included in Table 1, are neg- 
ligible compared to uncertainties resulting from continuum 
placement errors or line blending, which we discuss further 
below. The ultimate reason for the small Poisson uncer- 
tainties in these two components is that the S/N is high 
and the lines are relatively broad, spanning many pixels. 

Our derived Cf values are subject to a number of ad- 
ditional uncertainties. One such uncertainty results from 
line blending. If two or more components are blended to- 
gether, a direct estimation of the Cf value by the pixel- 
by-pixel method [eqn (2)] is not reliable. On the other 
hand, the results from the automatic fitting method are 
still reliable, especially for strong doublets, as long as they 
are not affected by self-blending. To evaluate the sys- 
tematic uncertainty in the latter method, we performed 
a line-deblcnding test in which we synthesized three ar- 
tificial C iv doublets [log(iV/cur 2 ) = 13.5, 14.0, 14.5; 
b = 35 km s _1 ; Cf = 0.5] separated by 100 km s _1 from 
each other, introduced them into an artificial spectrum 
with a S/N similar to what we obtain in practice, and de- 
termined their line parameters with MINFIT. We repeated 
this test 500 times. The distribution of resulting cover- 
ing factors for the stronger lines is very narrow and cen- 
tered very close to the input values; Cf = 0.51±0.06 and 
0.51±0.02 for doublets with log(iV/cm- 2 ) = 14.0 and 14.5. 
However, the spread in the measured value of Cf is very 
large for weak doublets with log(-/V/cm~ 2 ) = 13.5. The 



outcome does not change significantly if we change the in- 
put Cf values. Thus, the result from the automatic fitting 
method is still reliable, at least for strong doublets with 
\og(N/cm~ 2 ) > 14.0, even if there is line blending (but not 
self-blending). However, one must still attach a systematic 
error of < 10% to C f . 

Continuum fitting also introduces an uncertainty in the 
derived values of Cf because it affects the residual flux 
in the absorbed spectral regions. We estimated an upper 
limit to the uncertainty in C / resulting from continuum 
fitting errors through the following simulation: we shifted 
the normalized continuum level around the red member of 
the C iv doublet by the average noise level of the observed 
spectrum [cr(/) « 0.015] while keeping the continuum level 
around the blue doublet member fixed and we investigated 
what the new value of Cf turns out to be, under the as- 
sumption that Cf was originally equal to 1. From this 
exercise we found that we obtain Cf values 0.98, 0.93, 
0.82, 0.56, and 0.13 (instead of 1) for normalized residual 
fluxes at the center of blue doublet member of 0.1, 0.3, 
0.5, 0.7, and 0.9, respectively. Thus, we conclude that 
the Cf values for strong doublets (those with a normal- 
ized residual flux of < 0.5) are uncertain by <, 20%. We 
regard this as an upper limit to the uncertainty because 
the effect of shifting the continuum on the value of Cf is 
smaller if (a) we perturb the continuum around both dou- 
blet members, or (b) we assume that the original value of 
Cf is less than unity. In Table 3, we summarize upper and 
lower la errors on Cf values in the case that Cf is less 
than unity. These were computed by shifting the contin- 
uum level only around the red member. Column (1) gives 
the Cf value and columns (2) - (6) are the residual flux 
of the blue member of the doublet. We leave the entries 
blank for unphysical combinations of Cf and Rb. These 
tabulated values are the ones we adopt in our assessment 
of uncertainties. 

We also considered the uncertainty resulting from the 
convolution of the spectrum with the line spread function 
(LSF) of the spectrograph. This effect applies only to the 
pixel-by-pixel method (since automatic fitting convolves 
with the LSF) and was originally explored in the case of 
narrow doublets (widths comparable to the LSF width) by 
Ganguly et al. (1999). We performed the same analysis 
as Ganguly et al. but for broader doublets, because most 
doublets in our sample are much broader than the LSF. We 
found that for isolated lines that are much broader than 
the LSF, the Cf values are evaluated correctly within a 
window around the line center whose width is twice the 
Gaussian broadening parameter (26). This conclusion de- 
pends slightly on the line strength and width and is not 
valid if the doublet members arc blended with each other 
("self-blending"). Outside of this window, the Cf values 
are underestimated, and can reach in the far wings of 
the line. In these regions both the numerator and the de- 
nominator of equation (2) are very close to zero, therefore, 
the C f values obtained with the pixel-by-pixel method are 
sensitive to the pixel width, the noise, and the resampling 
and interpolation schemes used in the calculation. 

Considering all of the sources of error discussed above, 
the pixel-by-pixel method yields reliable results for our 
data in only very few cases. In the figures presented below, 
we plot the values of Cf from the pixel- by-pixel method 



6 



Misawa et al. 



only in the rare cases where their uncertainties are rea- 
sonably small. Otherwise, we rely solely on the value of 
Cf derived by fitting the line profiles (either manually or 
using minfit). 

4.2. Time Variability Analysis 

The detection of time variability of absorption line prop- 
erties such as strengths, profiles, and shift velocities is an- 
other reliable method for identifying intrinsic NALs. Vari- 
ability of absorption lines has been reported for a few 
Seyfert galaxies (e.g., Kraemer et al. 2002) and quasars 
(e.g., Gallagher et al. 2004; Hamann et al. 1997a). Re- 
cently, larger surveys looking for time variable NALs have 
been carried out at z ~ 2 (Narayanan et al. 2004) and 
at z < 1.5 (Wise et al. 2004), and they detected variable 
NALs in 2 of 8 (25%) and 4 of 15 (27%) quasars respec- 
tively. These fractions, however, are lower limits, since 
these studies used low- to medium-resolution spectra and 
are thus sensitive only to large variations. 

To search for spectral variability in HS 1603+3820, we 
overplotted selected regions of the spectra from the two 
epochs and inspected them by eye (cf, Hamann et al. 
1997a). This simple method is quite adequate because 
our spectra have a fairly high resolution and S/N so that 
even small changes are obvious to the eye. In summary, 
we found substantial variations in System A (Figure 4) 
but not in any other absorption system (Figure 5). We 
describe these changes along with other results in the fol- 
lowing section. Possible causes of variability are: (a) a 
change of ionization state, and (b) a motion of the absorb- 
ing gas across the line of sight. We discuss these possibil- 
ities further in §6. 




5300 5320 5340 

Wavelength (A) 

Fig. 4. — Normalized spectrum around C IV doublets in System A, 
in the 1st spectrum (thin histogram) and in the second spectrum 
(thick histogram). There is an obvious change in the absorption-line 
strength between the two epochs. Open and filled stars denote the 
covering factors for the first and second epoch respectively (from 
manual fitting) with lc errors based on the uncertainty from the 
line blending only (see discussion of additional uncertainties in §4.1 
and §5 of the text; Poisson errors are negligible in comparison). 

5. RESULTS: NOTES ON 8 C IV SYSTEMS 

We found 47 C iv, 14 Si iv, and 6 N v doublets, as well 
as another 54 single metal lines, such as Al II, Fe II, and 
Si II, between z a b s — 1-88 {v s h = 60,400 km s _1 ) and z a b s 



= 2.55 (v s h = —950 km s _1 ). For simplicity, we grouped 
the 121 components into 8 absorption systems, which are 
discussed below. For the first four systems, which have al- 
ready been discussed in M03, we performed both covering 
factor analysis and time variability analysis. For the other 
systems, we applied only covering factor analysis. 

System A (z abs = 2.4185 - 2.4455; v sh = 8,300-10,600 
km s -1 ) This system has many broad C IV doublets 
and some of them are self-blended, which makes the 
line profile very difficult to decompose into kinematic 
components. In this case, the results from both the 
pixel-by-pixel and the automatic fitting methods are 
unreliable. Therefore, we concentrate on the wave- 
length region at 5290-5315 A which can be decom- 
posed relatively simply into two kinematic compo- 
nents of C iv plus several narrow Si n lines, and 
employ the manual fitting method described in §4.1. 
If we assume full coverage, the best-fit model for the 
red (or blue) doublet member cannot reproduce the 
observed spectrum around the blue (or red) doublet 
member, as we illustrate in the top panel of Figure 2 
using the 2nd epoch spectrum as an example. Thus, 
we resort to models with partial coverage and we 
find coverage factors of 0.30 ± 0.02 and 0.45 ± 0.02 
for the first and second epoch spectra respectively 
(all sources of error included), i.e., we find a 5.3cr 
change in Cf. The uncertainties here are dominated 
by errors from line blending as estimated by the line 
blending analysis in § 4.1; uncertainties due to con- 
tinuum fitting (since Cf — 0.5 and Rb « 0.5; see 
Table 3) and Poisson noise are negligible in compar- 
ison (see §4.1 and Table 1). 

In the above analysis we assumed the same cover- 
ing factors for the two components. We were led 
to make this assumption for the following reason: 
we first fitted the velocity region between —300 and 
+250 km s _1 with a broad C iv component only and 
found that the required covering factor was the same 
as that obtained for the narrow component after fit- 
ting a 2-component (broad+narrow) model. After 
obtaining the best fit with the composite model we 
went back and relaxed the assumption that the two 
components had the same covering factors. We ex- 
plored the consequences by perturbing the covering 
factors of the two components about the best fit and 
adjusting the column density to compensate. From 
this exercise we found that, if the covering factors are 
not forced to be the same, they can deviate from the 
best fit by ±0.015 in the first epoch spectrum and 
by ±0.03 in the second epoch spectrum. If we treat 
these deviations as an additional error on Cf, we 
find that the change in Cf between the two epochs 
is 3.4cr, which is still fairly significant. 

To illustrate the change in Cf we plot the observed 
line profiles from the two epochs in the second and 
third panels on Figure 2 with the best fitting models 
for a fixed value of C/=0.45 superposed (all other 
parameters are free) . It is clear from this figure that 
the first-epoch profile cannot be fitted adequately us- 
ing the same value of Cf as the second-epoch profile. 
We also plot the best fitting model to the first-epoch 
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spectrum with Cf = 0.3 in the bottom panel. 

In Figure 3 we present a fit to the entire System A, 
including the heavily blended region around 5315- 
5345 A, for which we used the manual fitting method 
and assumed a fixed value of Cf = 1 . The parameters 
obtained from this fit are listed in Table 1, although 
they should be regarded as tentative because of their 
large uncertainties. Deriving the uncertainties for 
these model parameters would have been extremely 
tedious and is not central to the point of this paper. 

We identified 14 C iv doublets in this system, while 
D99 and M03 found 5 and 7 doublets, respectively. 
Lya, N v, and Si iv lines are also tentatively detected 
in our new spectrum. However, we cannot fit those 
profiles with minfit, because they are very broad 
and strong. The continuum fits in these regions are 
not reliable. We cannot even apply the interpola- 
tion technique described in §2, because echelle orders 
next to a given order are also affected by Lya forest. 

For comparison of line profiles between the 1st and 
2nd spectra, we also superimpose the two spectra 
in Figure 4. The strength of all doublets, including 
three narrow doublets, changes dramatically within 
1.28 years (0.36 years in the quasar rest frame). We 
also found that both the shift velocity and the line 
width of the broad component has changed, while 
such a change was not noted for component 1 . More 
specifically, the shift velocity of the broad compo- 
nent changed by ~ 210 km s _1 (the uncertainty in 
the shift velocity is ±5 km s^ 1 ), while the b param- 
eter changed by ~ 170 km s _1 (the uncertainty in 
b is ±8 km s _1 ). Unfortunately, we cannot study 



these parameters for the other components because 
of heavy line blending. The Si II A 1527 lines in 
System B at 5306-5314 A do not change between 
the two observations, which suggests our interpola- 
tion fitting method, described in §2 and in M03, is 
effective. Although M03 suggested that the narrow 
C iv doublets could be independent of the other lines 
in System A based only on the result of the cov- 
ering factor analysis, we find here that all 14 dou- 
blets change synchronously, and thus they are prob- 
ably produced in the same intrinsic system. These 
C iv doublets show a sudden increase of their total 
equivalent width (W obs « 10.4 A — » 19.1 A), which 
could mean that, in a later stage of evolution, many 
narrower components may blend with each other to 
make what resembles a BAL profile (Morris et al. 
1986; Richards et al. 2002; M03). Weymann et al. 
(1991) defined BALs as systems with continuous ab- 
sorption, spanning more than 2,000 km s _1 , within 
which at least 10% of flux must be absorbed. The 
velocity width of System A, Av ~ 2,300 km s" 1 , sat- 
isfies this definition, although the normalized resid- 
ual flux is more than 0.9 in some regions at the time 
of observation. Some examples of BALs are known 
to have different covering factors at different regions 
in the same system (Arav et al. 1999; Hutsemekers, 
Hall, & Brinkmann 2004). BAL systems also some- 
times have multiple-components, which are analo- 
gous to the components of System A (Korista et al. 
1993; Arav et al. 1999). These results support the 
hypothesis that the NAL cluster in System A could 
become a BAL, although additional observations will 
be necessary to confirm this. In §6.3 we use simula- 
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Fig. 5. — Same as Figure 4, but for System B, C, and D (binned by 10 pixels). The spectra from the two epochs are overplotted to illustrate 
the absence of variability in these three systems. 
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Fig. 6. — Velocity plots of various transitions in System B at z aba = 2.4785 (v abi = 5,400 km s _1 ). Vertical dotted lines with ID numbers 
are positions of absorption lines. Smooth lines are model fits by minfit. Filled stars with error bars are covering factors determined by 
minfit. Small dots with error bars denote covering factors evaluated by the pixel-by-pixcl method. In both cases, the error bars represent 
uncertainties from Poisson noise only; additional uncertainties are discussed in §4.1 of the text. 



tions to illustrate the plausibility of the suggestion 
that this system could be evolving into a BAL. 



System B (z abs = 2.4785; v sh = 5,400 km s" 1 ) This 

strong C iv system consists of 15 narrow doublets 
with b < 26 km s _1 . Some of them can be explained 
by pure thermal broadening (i.e., b < 7.4 km s _1 
at T gas = 40,000 K). We present the velocity plots 
of 12 metal line transitions, including the C iv and 
N v doublets, in Figure 6. In cases that MINFIT gives 
unphysical Cf values (i.e., Cf < or > 1), which 
happens for some components, probably because of 
heavy line blending, we assume Cf — 1 to estimate 
fit parameters. The covering factors of all doublets 
except for the component at z abs = 2.4804 (ID 18) 
are consistent with full coverage as determined both 
by profile fitting and the pixcl-by-pixel method. For 
the latter method, we do not plot Cf values if their 
error values are larger than 1.0. Component 18 of 
C iv A 1551 is likely to be influenced by the unidenti- 
fied blend apparent between components 17 and 18, 
thus the low Cf measured for component 18 of C iv 
A1551 is discounted. Component 5 of N v A1243 
may also be affected by a blend, since we were un- 
able to obtain a physical solution in this region for 
any reasonable continuum fit. 
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Fig. 7. — Same as Figure 6, but for System C at 
z abs = 2.5370 (v sh = 420 km s" 1 ). 

System C (z abs = 2.5370; v sh = 420 km s" 1 ) This is 
one of two systems within 1,000 km s _1 of the 
quasar. Besides the C iv doublet, there arc no 
other metal lines detected in the 2nd spectrum. We 
find neither partial coverage (Figure 7) nor time 
variability (see Figure 5) in spite of the very small 
velocity shift from the quasar. 

System D (z abs = 2.5532; v sh = -950 km s" 1 ) This is 
another system within 1,000 km s _1 of the quasar, 
but rcdshifted relative to the quasar. This system 
was resolved into 4 narrow doublets in M03. We are 
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able to analyze only three doublets in the 2nd spec- 
trum, because the bluest component of C iv A1548 
is affected by a detector defect. Although Si iv, as 
well as C IV, is identified in the system, neither shows 
partial coverage, if we consider the error in the con- 
tinuum fit (Figure 8). The Cf value of C iv compo- 
nent f in the 1st spectrum is small, but the upper 
la error would approach 1.0 once we consider con- 
tinuum level error, a(Cf) = 0.5 from Table 3. We 
also do not find time variability of the C iv doublets 
(see Figure 5). These results suggest that the ab- 
sorber covers the background source entirely in the 
direction of this line of sight. Nevertheless, we still 
cannot rule out an intrinsic origin for Systems C and 
D in the vicinity of the quasar (within 1,000 km s _1 ), 
especially because a number of intrinsic NALs have 
been detected at small velocities from quasars (e.g., 
Ganguly et al. 1999; Narayanan et al. 2004; Wise et 
al. 2004). We discuss these systems in §6. 
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Fig. 8. — Same as Figure 6, but for System D at 
z abs = 2.5532 (v sh = -950 km s" 1 ). 



System E (z abs = 1.8875; v sh = 60,400 km s" 1 ) Two 
C iv doublets with line widths of b = 16 and 27 
km s _1 are identified in the lowest redshift system 
in this study, which corresponds to v s h ~ 60,500 
km s _1 from the quasar (Figure 9). One of the two 
doublets (component 2) has a small covering factor, 
Cf — 0.5±0.2, and the errors from continuum level 
uncertainty (cr(Cf) — 0.1) and line blending (cr(C/) 
~ 0.06) are also very small (in § 4.1). Despite its 
large shift velocity, this system could be classified 
as an intrinsic system. To confirm this possibil- 
ity, additional evidence such as time variability is 
necessary. 
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Fig. 9. — Same as Figure 6, but for System E at 
Zabs = 1.8875 (v sh = 60,400 km s" 1 ). 



System F (z abs = 1.9644; v sh = 52,800 km s" 1 ) Out of 
eight C iv doublets in this system, five show full cov- 
erage (Figure 10). The other three doublets (com- 
ponents 1, 4, and 7) have small Cf values. One of 
them (component 4) remain a candidate for an in- 
trinsic doublet, even if we consider both the error 
sources from continuum level uncertainty (a(Cf) < 
0.1) and line blending (<r(C/) < 0.02) (see discussion 
in §4.1). 





1 2 3 


4 


;l 

5 f 


7 


1 n 






y 






: C IV 
: A154B 




1- 






1 '■ 






4 


5 £ 


7 




: C IV 
: A1551 

h 


y 


1- 









-200 



200 



velocity (km s 



Fig. 

z abs 



1 1 ! — Same as Figure 6, but for System F at 
1.9644 (v sh = 52, 800 km s" 1 ). 



System G (z abs = 2.0701; v sh = 42,600 km s" 1 ) This 
system must also cover the quasar entirely, because 
the line centers of the C iv A1548 transition are 
black (zero flux) which cannot happen in the case of 
partial coverage (Figure 11). The covering factors 
of components 1 and 2 are consistent with unity. 
Corresponding Si iv absorption is apparently de- 
tected at Av = —30 and km s _1 (see Figure 
11), though the components at Av > 50 km s" 1 
are N v doublets from System B. The components 
at Aw < —50 km s _1 for Si iv are also likely to 
arise from another system, because no correspond- 
ing C iv doublets are detected. Unfortunately, we 
cannot evaluate the covering factor of the Si IV dou- 
blet, because the Si iv A1394 line is blended with 
the Lya forest. 
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Fig. 11. — Same as Figure 6, but for System G at 
z ab3 = 2.0701 (v 3h = 42,600 km s" 1 ). 

em H (z abs = 2.2653; v sh = 24,300 km s" 1 ) A 
simple C iv system with two components, which 
was not detected in D99, is identified along with 
Si in A1207 (see Figure 12). The line center of Si in 
A1207 is well aligned with C iv doublet. Although 
the covering factor determined by minfit is smaller 
than unity (C/ = 0.5 ± 0.3) for component 1, we 
cannot rule out full coverage once we consider the 
uncertainties associated with setting the continuum 
level, <j(Cf) — 0.2, and with line blending (see dis- 
cussion in §4.1). This component, whose line width 
is comparable to the LSF, is affected by the con- 
volution error. Therefore the Cf values from the 
pixel- by-pixel method arc not reliable. There is a 
large discrepancy between the model fit and the data 
in component 2 of the C iv A1551 line. This suggests 
that the component 2 of the C iv A 1548 line could 
have been misidentified as C iv or could be affected 
by continuum fitting uncertainty. 
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Fig. 12. — Same as Figure 6, but for System H at 
z abs = 2.2653 (v sh = 24,300 km s" 1 ). 



6. DISCUSSION: ORIGIN OF TIME VARIABILITY 



Our main observational result is that System A shows a 
change of line strength within 0.36 years in the quasar rest- 
frame. No other system appears to be variable. The fact 
that System A is variable is compatible with the fact that 
it shows the signature of partial coverage and strengthens 
the case that the absorber is intrinsic to the quasar central 
engine. 

There are at least two possible origins of time variability, 

(a) a change of the ionization state in the gas clouds, and 

(b) motion of the absorbers across the line of sight to the 
background source(s). We discuss each of these possibili- 
ties below. Neither situation could arise if the absorbers 
are intervening clouds, unless they have high densities or 
sharp edges (Narayanan et al. 2004). 

6.1. Change of the Ionization State 

Supposing that the time variability is caused by a change 
in ionization state, we can place constraints on the elec- 
tron density and the distance from the continuum source, 
following Hamann et al. (1997b) and Narayanan et al. 
(2004). These limits should be regarded with some caution 
because of the assumptions involved in obtaining them. 
We outline these assumptions below. Changes in the ion- 
ization state of the absorber can be caused by (a) a de- 
crease in the ionizing flux causing recombination from the 
observed ionization state to the next lower state (C iv 
— > C in) or from the next higher state to the observed 
state (C v -t C iv), or (b) increase in the ionizing flux 
and subsequent photoionizations (C in — > C iv and C iv 
— ► C v). Unfortunately, with the available data we can- 
not assess whether the ionization state of the absorber 
of System A has changed or not. In order to do so we 
would have to resort to photoionization models, which 
would need to be constrained by absorption lines from 
a wide range of ionic species, which arc not included in 
our spectra. Moreover, we cannot detect any changes in 
the observed continuum because the flux scale of our spec- 
tra is not calibrated. Therefore, we adopt the simplifying 
assumptions of Hamann et al. (1997b) and Narayanan 
et al. (2004), namely: (1) the gas is very close to ion- 
ization equilibrium, (2) changes in the ionizing flux are 
small, and (3) C iv is dominant ionization state of C. 
Under these assumptions we can use the recombination 
time scale as a means of estimating the electron density 
[cqn. (1) in Hamann et al. 1997b; eqn. (1) in Narayanan 
et al. 2004; the recombination coefficient corresponds to 
a nominal gas temperature of 20,000 K, after Hamann et 
al. (1995)]. Using the rest-frame variability time scale 
{tvar = 0.36 years) as an upper limit to the recombina- 
tion time, we estimate the minimum electron density to 
be n e > 3 x 10 4 cm~ 3 . We can also evaluate the maximum 
distance of the absorber from the continuum source by 
adopting the continuum shape of Narayanan et al. (2004), 
and using their equation (3) to connect the ionization pa- 
rameter to the bolomctric luminosity [the bolometric lu- 
minosity of HS 1603+3820, based on the prescription of 
Narayanan et al. (2004) is L bo i = 2.5 x 10 48 ergs" 1 ]. 
Finally, we must also assume an ionization parameter of 
0.02, which is the optimal value for C iv to be the dom- 
inant ionization state of C (Hamann et al. 1995; 1997a). 
This yields r < 6 kpc. These limits on the electron density 
and the distance from the central source are very simi- 
lar to the results for C iv and N v systems seen in the 



Time- Variable Absorption Lines in HS 1603+3820 



11 



four quasars (n e > 3,400-40,000 cm -3 and r < 1.8-7 kpc) 
studied by Narayanan et al. (2004). 

In conclusion, we note that regardless of changes in ion- 
ization state, however, we can conclude with reasonable 
confidence that the absorber is moving across the line 
of sight because the covering factor of the background 
source(s) is changing. In the next section we consider the 
implications of the variable covering factor in detail. 

6.2. Motion of the Absorber 

In the case of motion of the absorbers, the variability 
time scale could correspond to (i) the crossing time of 
clumpy gas across the background source or (ii) the time 
scale on which the internal motion of the gas cloud changes 
the column density along the line of sight (e.g., Hamann 
ct al. 1997a; Wise et al. 2004). 

The variability time scale of System A is very short (~ 
130 days), suggesting that the observed changes could not 
be produced by internal motions of the gas cloud, as we 
will illustrate later in this section. Some intrinsic NALs 
detected in Narayanan et al. (2004) and Hamann et al. 
(1997a) also show very short variability time scales (t var 
= 0.28-0.40 years) in the quasar rest frame. In these sys- 
tems, the clumpy clouds could be moving across the line 
of sight to the continuum source. If this is the case for 
System A, we can derive some constraints on the size, 
crossing velocity, and distance of the absorber from the 
UV-emitting source using a simplified model as follows. 

At first, we assume that both the continuum source and 
BLR are background sources (e.g., Ganguly et al. 1999). 
In this case, we can estimate the size of the BLR to be 
Rblr ~ 3 pc by using the empirical relation between BLR 
size and the quasar luminosity that was found through 
reverberation mapping analysis [eqn (A3) in Vestergaard 
2002]. To apply this relation, we used the monochro- 
matic continuum luminosity derived from the B magni- 
tude of HS 1603+3820, measured by D99. This luminosity 
is strictly an upper limit because at the redshift of the 
quasar the B filter includes a contribution from the Lya 
and N V emission lines; nevertheless, the luminosity is 
overestimated only by a factor of a few. 

Supposing that the absorbers have sharp edges, we es- 
timate the crossing velocities as follows: 



>[C / (2)-C / (l)]x^ 



(3) 



where w cross is the crossing velocity perpendicular to our 
line of sight (in the context of a disk-wind model this is par- 
allel to the surface of the accretion disk), C/(l) and C/(2) 
are covering factors in the 1st and 2nd spectra, and t var 
is the time interval between the two observations. Since 
t var is really an upper limit to the variability time scale, 
equation (3) gives a lower limit to the crossing velocity. 
The crossing velocities of the two C iv components in 
System A (components 1 and 3 in the 1st spectrum, or 
components 1 and 2 in the 2nd spectrum), for which we 
can evaluate reliable C/ values in both spectra, would be 
Across > 10 6 km s _1 . This model is implausible, because 
the estimated speeds are much larger than the speed of 
light. Thus, we must conclude that the absorber covers 
primarily the UV continuum source and only a small part 
of the BLR. In fact, since the intensity of the broad C iv 
emission line at the location of System A is only about 



10% of the continuum, the absorber giving rise to Sys- 
tem A need not cover the BLR at all. 

Therefore, we estimate the velocities by assuming that 
only the continuum source lies behind the absorbing clouds 
(e.g., Wise et al. 2004). We take 5 R s = 10 GM BH /c 2 
(where Rs and Mbh are the Schwarzschild radius and the 
mass of black hole) as the continuum source size within 
which most of the UV continuum radiation is expected 
to originate. The central black hole mass is estimated by 
using an empirical relation that connects it to the UV lu- 
minosity and the width of the C IV emission line [eqn (8) in 
Vestergaard (2002)]. The UV luminosity was estimated as 
described above, and the FWHM of the C IV emission line 
(7,500 km s _1 ) was measured from the spectrum of D99. 
The resulting black hole mass is 4 x 10 10 M©, but we note 
that this value should be regarded with some caution be- 
cause (a) the scatter about the empirical relation we have 
used can be as high as an order of magnitude, and (b) the 
luminosity we have adopted is strictly an upper limit to 
the true value. Therefore, we believe that the black hole 
mass can be as low as a few xlO 9 M . This model pro- 
duces reasonable values of the crossing velocities, namely 
Vcross > 8,000 km s^ 1 , using the higher black hole mass. 
We may also use the black hole mass to estimate an upper 
limit to the distance of the absorber from the background 
source, by requiring that the observed shift velocity does 
not exceed the escape velocity at that radius. From this 
consideration we obtain r < 3 pc. 

The above estimate of the crossing velocity can be used 
to argue against internal motions of the absorber as the 
origin of the change in the covering factor. Assuming that 
the speed of internal motions is comparable to the speed of 

sound, c s s=s 10 T^ 2 km s _1 (where T4 is the temperature 
in units of 10 4 K), the gas must be very hot to allow such 
fast internal motions. In particular, by setting v cross = c s 
we obtain T gas ~ 6 x 10 8 K, which is incompatible with 
the observed ionic species. We note, in conclusion, that 
using the size of the UV continuum source from the pre- 
vious paragraph and the speed of sound for a 10 4 K gas, 
a typical variability time scale due to internal gas motions 
would be of order 100 years. 

We can derive more stringent, albeit more specula- 
tive, constraints on the location of the absorber, if we 
adopt the dynamical model of Murray et al. (1995) for 
the radial and azimuthal motion of parcels of gas in an 
accretion-disk wind. Following this model, we assume 
that the velocity of the gas in the radial direction follows 
v r {r) ~ Uoo(r/)(l — Tj/r), where 77 is the launch radius 
of a specific gas parcel and t , 00 (?'/) is its terminal speed. 
Combining this with conservation of angular momentum 
and the expression for the terminal speed given in Chiang 
& Murray (1996) we obtain the following relation for the 
radius of the gas parcel relative to its launch radius: 

~v r (r) ' 



— w 0.2 

r f 



,(r) 



1 



(4) 



where v v {r) is the azimuthal speed of the parcel. To con- 
nect this expression to the observations, we identify the 
azimuthal velocity with the crossing velocity of the contin- 
uum source inferred from observations, i.e. v v (r) — v cross . 
Moreover, we can take the velocity of System A to be the 
projection of the radial velocity along the line of sight, i.e., 
v r (r) sin i — vs ys A ■ If the inclination of the accretion disk 



12 



Misawa et al. 



is not very small (sini > 0.2), we can use equation (4) 
to obtain an upper limit on r/rf. This limit will depend 
directly on the adopted mass of the black hole, since that 
determines the inferred crossing velocity. In view of the 
caveats associated with the estimated black hole mass, we 
choose to use the lowest plausible value, namely 10 9 M Q , 
which yields the most generous upper limit on the radius of 
the absorber, namely r/rf < 10. The launch radius of the 
wind is expected to be within the inner, radiation-pressure 
supported part of the accretion disk, which for our inferred 
black hole mass and assuming an accretion rate close to the 
Eddington limit, is 5 x 10 16 cm (Shakura & Sunyaev 1973; 
this is rather insensitive to the black hole mass) . Thus the 
constraint on the radius translates to r < 0.2 pc. We also 
note that the conclusion that r/rf < 10 is in good agree- 
ment with the simulations of Proga et al. (2000), which 
show the wind fragmenting into filaments fairly close to its 
launch radius. 7 



6.3. Evolution to a BAL 

The total equivalent width of System A has doubled 
within 0.36 years in the quasar rest frame. This surprising 
growth is probably caused by (i) increase of column densi- 
ties, and/or (ii) increase of the covering factors of the C iv 
doublets in the system, although the Doppler parameters 
could also be partially responsible. If the equivalent width 
continues to increase steadily, System A would satisfy all 
the requirements for a BAL classification. 

In order to predict future profiles, we synthesize spectra 
by changing the column densities or the covering factors of 
the C iv doublets in the system. Unfortunately, we cannot 
evaluate the column density and covering factor of each 
doublet exactly in the observed spectrum, since the C iv 
doublets are heavily blended with each other or with Si n 
lines from System B. Therefore, we convert the column 
densities of the C iv doublets in the 2nd spectrum in Table 
1 (evaluated assuming Cf = 0.45 or 1.0) to values assum- 
ing Cf = 0.7 (corresponding to the depth of the deepest 
C iv component; we adopt this value for all components 
for simplicity), and use them below. We present simulated 
line profiles which are synthesized by multiplying the col- 
umn densities by 0.5, 1.0, 2.0, and 5.0 (Figure 13; here, 
covering factors are fixed), or by changing the covering fac- 
tors, Cf — 0.3, 0.5, 0.7, and 0.9 (Figure 14; here, column 
densities are fixed). Thick lines arc the observed profiles, 
while thin lines are synthesized models. For three mod- 
els (i.e., N 2 /N 1 = 2.0, 5.0, or C f = 0.9), the continuum 
flux is absorbed by at least 10% over a velocity range of > 
2,000 km s _1 , which satisfies the definition of a BAL. If the 
Doppler parameters increase simultaneously, the blended 
C iv components make a smoother profile, whose appear- 
ance would be similar to those of BALs. With additional 
spectra taken at regular intervals, we could probe the ori- 
gin of time variability further (i.e., distinguish between a 
change of column density and a change in covering factor) 
by comparing the observed line profiles with our simula- 
tions. 

7 If we had chosen the higher black hole mass of 4 X 10 10 Mq, we would 



T 



C, = 0.7 (const.) 




5300 5320 5340 

Wavelength (A) 

Fig. 13. — Model profiles of System A, synthesized by multiplying 
column densities of the C IV components by factors of N2/N1 = 0.5, 
1.0, 2.0, and 5.0 (top to bottom). The covering factors are assumed 
to be constant (Cf = 0.7 in all cases). The original column densities 
were estimated, assuming that Cf = 0.7 for all C IV doublets in the 
observed spectrum (except for the Si II lines in System B). The thick 
line corresponds to the original (observed) profile, while thin lines 
are model profiles for different column densities. The total velocity 
width of the system is about 2,300 km s _1 . Therefore, the model 
profile is classified as a BAL, if its normalized flux is weaker than 0.9 
(horizontal dotted line) throughout the system. Two of the models 
(JV2/JV1 = 2.0 and 5.0) satisfy all the BAL criteria. 



T 



log N = const. 




5300 5320 5340 

Wavelength (A) 

Fig. 14. — Same as Figure 13, but synthesized by changing cov- 
ering factors, Cf = 0.3, 0.5, 0.7, and 0.9 (top to bottom). Column 
densities are fixed. One of these models (Cf = 0.9) satisfies the 
criteria for a BAL. 

6.4. Physical States of System C and D 

In Systems C and D, we detect neither time variability 
nor partial coverage with our spectra taken over a 1.28- 
year interval of observation, even though these would ap- 
pear to be promising candidates for intrinsic systems be- 
cause of their small shift velocities. 

There are at least three candidates for the origin of these 
systems. They can be produced in either (i) the quasar 
host galaxy, (ii) intervening galaxies in the vicinity of the 
quasar that are members of the same cluster (or group) 
of galaxies, or (iii) gas clouds intrinsically associated with 

have obtained r/rf ~ 1 and r < 0.02 pc. 
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the quasar that appear to be ejected (System C) or in- 
falling (System D) ; these gas parcels can be part of a high 
altitude stream in an accretion disk wind (which may be 
circulating) as simulated by Proga et al. (2000). 

However, since we do not see the hallmarks of intrin- 
sic absorption in cither of these two systems (e.g., partial 
coverage or time variability), we cannot prove that they 
are intrinsic, nor can we rule this out. We find it unlikely 
that System D is made up of undisturbed gas in the host 
galaxy because of its large velocity shift. The fact that 
the absorption troughs in System D are black, suggests 
that this absorber covers both the continuum source and 
the broad-emission line region almost completely, i.e., its 
transverse size could be a few pc (cf, §6.2). This large size 
could also be responsible for the lack of variability, if this 
system is intrinsic. 

6.5. Other Potentially Intrinsic Systems 

Among the 8 C iv systems, four (Systems E - H) have 
been studied based on only covering factor analysis. In 
cases where the measured Cf value is small, the errors 
of Cf from continuum level uncertainty, as listed in Ta- 
ble 3, are small, and such components can be classified 
as intrinsic doublets: component 2 in System E (C/ = 
0.5±0.2) and component 4 in System F (Cf — 0.1 with 
o~(Cf) smaller than 0.1) (here errors from the line blending 
are negligible). These systems are potentially interesting 
because of their high velocities relative to the quasar: v s h 
~ 60,400 km s" 1 (System E) and 52,800 km s" 1 (System 
F). Therefore, it is very important to monitor them further 
for time variability. We already know of several examples 
of high- velocity systems, which are variable. These include 
the systems in Q0151+048 with v sh ~ 28,000 km s" 1 , 
in PG0935+417 with v sh ~ 52,000 km s" 1 (Narayanan 
et al. 2004), in Q2343+125 with v sh ~ 24,000 km s _1 
(Hamann et al. 1997a), and in PG2302+029 with v sh ~ 
56,000 km s" 1 (Jannuzi 2002). 

7. CONCLUSIONS AND FUTURE WORK 

We observed the quasar HS 1603+3820 for a second time 
with Subaru+HDS after an interval of 1.28 years from the 
previous observation, in order to perform time variability 
analysis for the C IV systems at v s h — —950 km s _1 - 
10,600 km s _1 from the quasar. We also applied cover- 
ing factor analysis using two methods, the pixel-by-pixcl 
method and Voigt profile fitting. Our main results are: 

(1) Only System A (v sh ~ 8,300-10,600 km s" 1 ), which 
was identified as an intrinsic system by covering fac- 
tor analysis in M03, shows strong time variability of 
the line strength and covering factor. 

(2) If a change of the ionization state causes the variabil- 
ity in System A, we can place constraints on the elec- 
tron density (n e > 3 x 10 4 cm~ 3 ) and the absorber's 
distance from the continuum source (r < 6 kpc). 



(3) If gas motion across the background UV source 
causes the variability in System A, the crossing ve- 
locity and the distance from the continuum source 
are estimated to be v cross > 8,000 km s _1 and 
r < 3 pc, respectively. In this case, the BLR cannot 
contribute significantly to the background UV flux 
(in fact, it need not contribute at all). If we adopt 
the dynamical model of Murray et al. (1995), we 
obtain an additional relation for the radius of the 
gas parcel relative to its launch radius, which gives 
a more strict, although model-dependent, constraint 
of r < 0.2 pc. 

(4) We argue that the appearance of System A would be 
fairly similar to that of a BAL system, if cither the 
column densities of the C IV lines increase by factors 
of > 2, or if their covering factors become larger than 
0.9. 

(5) Systems C and D (v sh ~ 420 km s" 1 and -950 
km s _1 ) do not show any time variability, in spite of 
their small velocity distance from the quasar. They 
could still be intrinsic systems but a more stringent 
test will have to await future observations 

(6) Two systems (System E and F; v s h — 60,400 and 
52,800 km s _1 ) show suggestive evidence of an in- 
trinsic origin based on only covering factor analysis. 
To confirm this, we would need other evidence for 
intrinsic NALs, such as time variability. 

In order to investigate the nature of System C further, 
there are at least two possible future observations. If we 
observe this quasar repeatedly at the appropriate time in- 
tervals (e.g., once a year), we may detect time variability 
in System C. This observation would have the additional 
benefit that we could detect the formation of a BAL in Sys- 
tem A as we have speculated. Since the equivalent width of 
System A has doubled in only ~ 1.3 years, a few more ob- 
servations may show dramatic evolution. We could also di- 
agnose variations in the ionization state of the absorber by 
comparing observed spectra with our simulations. Spec- 
tropolarimctric observations would also be worthwhile be- 
cause the degree of polarization in the troughs of intrin- 
sic NALs could be much larger than those in intervening 
NALs, indicating scattering around the absorber. We be- 
lieve that HS 1603+3820, one of the best laboratories for 
the investigation of intrinsic NALs, is an excellent target 
to monitor over a long period of time. 
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Table 1 

Metal lines detected in 1st and 2nd spectra 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

Line ID X obs z abs V b C } log N X obs z abs V b C f log AT 

(A) (kms- 1 ) (kms- 1 ) (cm- 2 ) (A) (kms" 1 ) (kms" 1 ) (cm- 2 ) 



C iv A1548: 

I 5292.6 

2 

3 5301.6 

4 5316.5 

5 5319.8 

6 5320.4 

7 5322.8 

8 

9 5326.1 

10 5327.6 

II 5327.6 

12 5329.4 

13 

14 5331.7 

15 5334.2 

Al nXWTT: 

1 5811.2 

2 5812.0 

3 5813.3 

4 5813.8 

5 5815.2 

6 5816.3 

7 5816.7 

C IVA1548: 

1 5381.8 

2 5382.2 

3 5382.8 

4 5383.0 

5 

6 

7 5383.8 

8 5383.9 

9 5384.5 

10 5384.5 

11 5385.0 

12 5385.4 

13 5385.6 

14 5386.1 

15 5386.3 

16 5387.1 

17 5387.4 

18 5388.4 

19 5388.6 

Fe IIA1608: 

1 5595.1 

Si HA1527: 

1 5302.9 

2 5303.7 

3 5304.6 

4 5306.4 

5 5307.4 

6 5310.2 

7 5310.8 

8 5312.5 



System A : z abs = 2.4302 



2.4185 

2.4246 
2.4340 
2.4361 
2.4365 
2.4381 

2.4402 
2.4412 
2.4412 
2.4423 

2.4438 
2.4454 



10646 

10106 
9287 
9104 
9069 
8929 

8746 
8659 
8659 
8563 

8433 
8293 



86 ± 2 a 

238 ± 5 a 
60.3° 
78. 5° 
27.1° 
50. 7° 

289. 3° 
9.9 b 
42.1° 
38. 9 b 

59A b 
54.8° 



0.31 ± 0.01 a 
0.31° 



14.59 ± 0.04 a 

14.59 ± 0.01" 
13.60° 
13. S2 b 
13.45° 
12. 84° 

13. 97 6 
13.23° 
13.47° 
13.29° 

13.34° 
13.54° 



5292.6 
5298.4 



2.4185 
2.4222 



10639 
10319 



88 ± 1" 
406 ± 8 a 



0.44 ± 0.01° 
0.44 a 



14.62 ± 0.01 a 
14.76 ± 0.01° 



System B 



Zabs 



2.4781 
2.4786 
2.4794 
2.4797 
2.4805 
2.4812 
2.4814 

2.4762 
2.4764 
2.4768 
2.4770 



2.4775 
2.4775 
2.4779 
2.4779 
2.4783 
2.4785 
2.4787 
2.4790 
2.4791 
2.4796 
2.4798 
2.4804 
2.4806 

2.4786 

2.4737 
2.4742 
2.4748 
2.4757 
2.4765 
2.4782 
2.4786 
2.4798 



5461 
5418 
5349 
5323 
5254 
5194 
5177 

5627 
5606 
5571 
5559 



5517 
5511 
5477 
5475 
5447 
5428 
5413 
5386 
5378 
5334 
5317 
5262 
5250 

5418 

5840 
5797 
5746 
5668 
5608 
5452 
5418 
5323 



11.9 ± 0.8 
9.3 ± 0.2 

4.8 ± 1.6 

5.7 ± 0.6 

2.9 ± 0.3 

3.8 ± 0.7 
6.2 ± 1.6 

8.5 ± 0.3 C 

6.6 ± 0.2 C 
23.2 ± 1.0° 
2.2 ± 1.3° 



16.5 ± 2.0 
7.2 ± 0.6 C 
14.9 ± 5.5 

4.1 ± 1.4 
13.5 ± 1.4° 

5.0 ± 3.0 
12.8 ± 6.2 
7.7 ± 1.7 

27.2 ± 0.5 C 

6.2 ± 1.2 
8.2 ± 0.3 
18.2 ± 2.4 
4.7 ± 0.2 

4.1 ± 0.5 

22.9 ± 4.8 C 
3.1 C 
4.8° 
0.1 C 
0.8° 
1.0° 



0.95 ± 0.03 

1.01 ± 0.15 
0.85 ± 0.56 

0.99 ± 0.47 
0.98 ± 0.01 
0.89 ± 0.07 

1.29 ± 1.50 
0.98 ± 0.02 
0.59 ± 0.77 
0.97 ± 0.02 



11.99 ± 0.02 
12.88 ± 0.02 
11.23 ± 0.08 
11.75 ± 0.03 
11.95 ± 0.03 
11.71 ± 0.05 
11.55 ± 0.08 

13.36 ± 0.01° 
14.08 ± 0.02 c 
13.49 ± 0.02 c 
12.27 ± 0.09° 



13.67 ± 0.18 
14.16 ± 0.04 c 
13.59 ± 0.09 
13.10 ± 0.43 
14.03 ± 0.05 c 
13.09 ± 0.70 
13.99 ± 0.19 
13.93 ± 0.06 
13.90 ± 0.01° 
12.73 ± 0.57 
13.75 ± 0.03 
13.25 ± 0.65 
13.97 ± 0.06 

13.69 ± 0.03 



21.8 
23.6 
11.2 
13.8 
14.2 



± 
± 
± 
± 
± 

11.2 ± 0.1' 

12.3 ± 0.4' 



13.06 
13.14 
12.95 
13.63 



0.07° 
0.02° 
0.04 c 
0.01° 



13.04 ± 0.03° 
13.24 ± 0.02 c 
14.07 ± 0.02 c 
13.17 ± 0.01° 



5316.8 


2.4342 


9269 


70.7° 










13.91° 


5319.9 


2.4362 


9095 


77.7° 










14.20° 


5320.5 


2.4366 


9060 


23.2° 










13.36° 


5322.8 


2.4381 


8929 


48.1° 










13.53° 


5324.1 


2.4389 


8859 


31.5° 










13.15° 


5325.8 


2.4400 


8764 


84.1° 










13.91° 


5327.6 


2.4412 


8659 


12.3° 










13.54° 


5327.6 


2.4412 


8659 


41.1° 










13.80° 


5329.4 


2.4423 


8563 


51.8° 










13.96° 


5330.7 


2.4432 


8485 


24.1° 










12.94° 


5331.8 


2.4439 


8424 


49.0° 










13.70° 


5334.3 


2.4455 


8285 


82.1° 










13.99° 


.4785 


5811.3 


2.4782 


5452 


10.6 ± 1.1 










11.88 ± 0.04 


5812.0 


2.4786 


5418 


9.3 ± 0.3 










12.81 ± 0.03 


5813.3 


2.4794 


5349 


1.3 ± 3.5 










11.03 ± 0.14 


5813.8 


2.4797 


5323 


4.8 ± 0.9 










11.62 ± 0.05 


5815.2 


2.4805 


5254 


2.8 ± 0.5 










11.91 ± 0.04 


5816.3 


2.4812 


5194 


3.1 ± 1.2 










11.62 ± 0.10 


5816.5 


2.4813 


5185 


5.7 ± 3.5 










11.36 ± 0.19 


5381.8 


2.4762 


5625 


8.6 ± 0.6 


1 


54 


± 


0.61 


13.07 ± 0.21 


5382.1 


2.4764 


5608 


7.3 ± 0.3 





98 


± 


0.02 


14.06 ± 0.05 


5382.8 


2.4768 


5573 


14.2 ± 2.0° 










13.31 ± 0.06° 


5383.1 


2.4770 


5556 


3.6 ± 2.1 c 










12.51 ± 0.22° 


5383.4 


2.4772 


5539 


9.7 ± 3.9° 










11.51 ± 9.99° 


5383.5 


2.4773 


5530 


25.0 ± 2.0° 










13.27 ± 0.23° 


5383.8 


2.4775 


5513 


9.0 ± 0.5 





99 


± 


0.02 


14.22 ± 0.05 


5384.5 


2.4779 


5478 


12.5 ± 0.9 


1 


00 


± 


0.07 


13.72 ± 0.06 


5385.1 


2.4783 


5444 


12.4 ± 1.4 





98 


± 


0.03 


14.02 ± 0.05 


5385.6 


2.4786 


5418 


12.9 ± 2.9 





85 


± 


0.08 


14.00 ± 0.09 


5386.2 


2.4790 


5383 


26.1 ± 2.1 





99 


± 


0.03 


14.20 ± 0.06 


5387.3 


2.4797 


5323 


16.4 ± 1.8 





69 


± 


0.29 


13.68 ± 0.26 


5387.4 


2.4798 


5314 


7.0 ± 0.9 


1 


02 


± 


0.19 


13.49 ± 0.16 


5388.3 


2.4804 


5263 


11.3 ± 4.4 





16 


± 


0.02 


14.05 ± 0.36 


5388.6 


2.4806 


5246 


4.7 ± 0.4 





93 


± 


0.03 


14.06 ± 0.15 


5595.1 


2.4786 


5418 


6.3 ± 0.5 










13.60 ± 0.03 


5302.9 


2.4737 


5840 


19.0 ± 6.2° 










13.04 ± 0.06° 


5303.7 


2.4742 


5797 


23.7 ± 3.8° 










13.10 ± 0.07° 


5304.6 


2.4748 


5746 


19.1 ± 4.0° 










12.95 ± 0.08° 


5306.4 


2.4757 


5668 


10.9 ± 0.8° 










13.61 ± 0.03° 


5307.6 


2.4765 


5599 


12.3 ± 1.1° 










13.03 ± 0.06° 


5310.0 


2.4782 


5461 


14.2 ± 2.3° 










13.24 ± 0.06° 


5310.8 


2.4786 


5418 


11.2 ± 0.3° 










14.07 ± 0.03° 


5312.5 
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5323 


12.3 ± 1.5° 
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Table 1 — Continued 



03 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

Line ID \ obs z abs V b Cj log N X obs z abs V b C; log N 

(A) (kms- 1 ) (kins" 1 ) (cm- 2 ) (A) (kms" 1 ) (kms" 1 ) (cm- 2 ) 



9 5313.7 2.4805 5254 4.7 ± 0.2 C ... 13.45 ± 0.01 c 5313.9 2.4805 5246 5.6 ± 0.2 C ... 13.41 ± 0.03° 

10 5314.8 2.4812 5194 7.2 ± 0.5 C ... 13.22 ± 0.03 c 5314.8 2.4812 5194 7.8 ± 0.9° ... 13.22 ± 0.05 c 

11 5315.1 2.4815 5177 11.7 ± 0.4 C ... 13.43 ± 0.02 c 5315.2 2.4815 5168 11.5 ± 0.4 C ... 13.49 ± 0.02° 



System C : z abs = 2.5370 

C IVA1548: 

1 5475.8 2.5369 432 11.6 ± 0.3 1.10 ± 0.17 13.30 ± 0.09 5475.8 2.5369 432 11.4 ± 0.5 0.88 ± 0.16 13.44 ± 0.11 



System D : z abs — 2.5532 

(J IVA1548: 

1 5497.9 2.5512 -778 4.8 ± 0.7 0.37 ± 0.21 13.07 ± 0.35 

2 5500.1 2.5526 -896 11.0 ± 0.6 0.91 ± 0.35 13.15 ± 0.20 5500.1 2.5526 -896 9.5 ± 1.4 0.59 ± 0.34 13.34 ± 0.35 

3 5500.9 2.5531 -939 13.1 ± 0.2 0.98 ± 0.01 14.28 ± 0.01 5501.0 2.5532 -947 14.4 ± 1.1 0.99 ± 0.03 14.17 ± 0.07 

4 5501.5 2.5535 -972 8.1 ± 0.2 0.99 ± 0.00 15.05 ± 0.08 5501.5 2.5535 -972 8.9 ± 1.0 0.97 ± 0.02 14.82 ± 0.25 



a Parameters are obtained by manual fitting method with Cf left as a free parameter. The uncertainties given here are the result of Poisson errors only 
and were evaluated by scanning parameter space, as discussed in §4.1 of the text. The covering factors for components 1 and 3 (in the 1st spectrum) and 1 
and 2 (in the 2nd spectrum) were required to be the same. 

b Paramctcrs arc obtained by manual fitting, assuming Cf — 1.0. The uncertainties are discussed in §4.1 of the text. 

c Parameters are obtained by automatic fitting, assuming Cf — 1.0. The errors quoted here reflect only the uncertainties due to Poisson noise. 



Time- Variable Absorption Lines in HS 1603+3820 



Table 2 

Metal lines in the 2nd spectrum 



(i) 

Line ID 



(2) 

(A) 



(3) 



(4) 
V 

(km s -1 ) 



(5) 
b 

(km s _1 ) 



(6) 



(7) 

log N 

(cm -2 ) 



System B : z abs = 2.4785 



IlAl335: 

1 4639.4 2.4764 5608 4.8 ± 0.9 

2 4641.6 2.4781 5461 13.0 ± 1.4 

3 4642.2 2.4785 5427 10.4 ± 3.3 

4 4642.4 2.4787 5409 5.4 ± 1.8 

5 4642.8 2.4790 5383 17.8 ± 9.7 

6 4643.5 2.4795 5340 6.7 ± 2.1 

7 4643.8 2.4797 5323 6.0 ± 0.6 

8 4644.8 2.4805 5254 6.0 ± 0.6 

9 4645.8 2.4812 5194 4.0 ± 1.8 

10 4645.9 2.4813 5185 4.4 ± 1.5 

11 4646.3 2.4816 5159 5.6 ± 2.3 

N VA1239: 

1 4308.0 2.4775 5513 7.9 ± 1.3 0.70 ± 0.49 

2 4308.5 2.4779 5478 7.7 ± 4.3° 

3 4308.9 2.4782 5452 13.5 ± 4.0° 

4 4309.2 2.4785 5427 7.7 ± 1.1° 

5 4309.5 2.4787 5409 ... b 

6 4309.9 2.4790 5383 9.0 ± 1.5" 

O IA1302: 

1 4529.6 2.4785 5427 8.3 ± 0.3 

2 4533.1 2.4812 5194 12.2 ± 5.4 

Si HA1260: 

1 4381.4 2.4761 5634 4.6 ± 1.5 

2 4381.7 2.4764 5608 3.1 ± 0.7 

3 4383.6 2.4779 5478 3.2 ± 1.4 

4 4383.9 2.4781 5461 10.7 ± 1.2 

5 4384.4 2.4785 5427 2.6 ± 8.1 

6 4384.6 2.4787 5409 10.2 ± 3.6 

7 4385.5 2.4794 5349 6.4 ± 1.6 

8 4385.9 2.4797 5323 5.3 ± 0.5 

9 4386.9 2.4805 5254 5.4 ± 0.4 

10 4387.8 2.4812 5194 4.9 ± 0.7 

11 4388.0 2.4814 5177 3.2 ± 0.9 

Si HIA1207: 

1 4194.9 2.4769 5564 28.4 ± 4.5 

2 4196.1 2.4779 5478 55.3 ±35.3 

3 4196.9 2.4786 5418 8.8 ± 8.9 

4 4197.4 2.4790 5383 12.4 ± 5.2 

5 4198.0 2.4795 5340 7.3 ± 3.4 

6 4198.3 2.4797 5323 32.7 ± 6.0 

7 4198.3 2.4797 5323 2.6 ± 0.4 

8 4199.2 2.4805 5254 7.4 ± 1.2 

9 4200.1 2.4812 5194 4.1 ± 1.4 

10 4200.3 2.4814 5177 5.2 ± 1.9 

Si IVA1394: 

1 4845.1 2.4763 5616 12.7 ± 5.7 

2 4845.2 2.4764 5608 4.5 ± 0.5 

3 4846.8 2.4775 5513 6.6 ± 0.8 

4 4847.3 2.4779 5478 3.3 ± 0.6 

5 4847.8 2.4782 5452 20.6 ± 3.2 

6 4848.3 2.4786 5418 9.6 ± 0.8 

7 4848.9 2.4790 5383 14.4 ± 1.4 

8 4849.7 2.4796 5332 24.8 ± 2.3 

9 4850.0 2.4798 5314 4.6 ± 0.9 

10 4851.0 2.4805 5254 5.8 ± 0.5 

Si IVA1394: 

1 4952.2 2.5531 -939 10.0 ± 1.1 0.77 ± 0.42 

2 4952.7 2.5535 -972 8.5 ± 0.4 0.94 ± 0.07 

(J IVA1548: 

1 4470.1 1.8873 60471 27.1 ± 

2 4470.9 1.8878 60422 16.1 ± 

(J IVA1548: 

1 4587.3 1.9630 52985 11.2 ± 2.3 0.47 ± 0.19 

2 4587.8 1.9633 52956 7.1 ± 0.7 0.96 ± 0.05 

3 4587.9 1.9634 52946 7.7 ± 0.5 0.98 ± 0.06 

4 4588.5 1.9638 52907 3.0 ± 2.6 0.11 ± 0.02 

5 4589.5 1.9644 52848 8.3 ± 0.5 0.96 ± 0.03 

6 4589.8 1.9646 52828 11.7 ± 3.3 0.69 ± 0.40 

7 4590.2 1.9649 52799 8.6 ± 1.1 0.76 ± 0.15 



System D : z abs = 2.5532 



13 


28 


± 


0.06 


14 


04 


± 


0.05 


14 


28 


± 


0.13 


13 


87 


± 


0.23 


13 


09 


± 


0.17 


12 


97 


± 


0.10 


13 


75 


± 


0.06 


13 


96 


± 


0.10 


14 


25 


± 


0.66 


13 


76 


± 


0.15 


12 


82 


± 


0.12 


13 


41 


± 


0.49 


12.48 ± 0.21° 


12.90 ± 0.11" 


12.69 ± 0.05° 

b 


12.65 ± 0.12° 


14 


34 


± 


0.03 


13 


43 


± 


0.15 


12 


03 


± 


0.08 


12 


49 


± 


0.08 


12 


18 


± 


0.11 


13 


04 


± 


0.03 


15 


96 


± 


0.49 


12 


90 


± 


0.31 


12 


11 


± 


0.08 


12 


82 


± 


0.05 


13 


02 


± 


0.06 


13 


09 


± 


0.10 


12 


41 


± 


0.07 


14 


34 


± 


0.38 


13 


79 


± 


0.24 


14 


05 


± 


2.42 


12 


75 


± 


0.30 


12 


95 


± 


0.28 


13 


10 


± 


0.14 


14 


50 


± 


0.68 


13 


24 


± 


0.24 


12 


48 


± 


0.11 


12 


34 


± 


0.11 


12 


10 


± 


0.34 


13 


10 


± 


0.05 


12 


43 


± 


0.04 


12 


66 


± 


0.05 


13 


01 


± 


0.05 


13 


08 


± 


0.04 


12 


89 


± 


0.03 


12 


99 


± 


0.04 


12 


74 


± 


0.07 


13 


17 


± 


0.06 



System E 



= 1.8875 



12.91 ± 0.30 
13.39 ± 0.07 



2.9 
3.2 



0.95 
0.50 



0.28 
0.23 



Systc 



13.80 ± 0.17 
13.70 ± 0.30 



F : z ab3 = 1.9644 



13.50 ± 0.26 
13.86 ± 0.07 
13.76 ± 0.06 
15.82 ± 3.57 
13.98 ± 0.06 
13.34 ± 0.33 
13.48 ± 0.13 



Misawa et al. 
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(1) (2) (3) (4) (5) (6) (7) 

Lino ID \ obs z abs V b C f log N 

(A) (km s- 1 ) (km s- 1 ) (cm- 2 ) 



8 4590.4 1.9650 52789 7.2 ± 1.5 0.73 ± 0.69 13.06 ± 0.51 



System G : z abs = 2.0701 

C IVA1548: 

1 4752.5 2.0697 42643 6.6 ± 1.1 1.29 ± 0.40 13.38 ± 0.23 

2 4753.0 2.0700 42614 19.4 ± 1.8 1.04 ± 0.05 14.16 ± 0.05 

Si IVA1403: 

1 4306.2 2.0698 42633 13.2 ± 2.5 ... 13.10 ± 0.06 

2 4306.6 2.0701 42604 6.8 ± 2.1 ... 12.75 ± 0.11 



System H : z ab3 = 2.2653 
IVA1548: 

1 5055.2 2.2652 24357 5.8 ± 1.8 0.49 ± 0.30 13.23 ± 0.33 

2 5055.5 2.2654 24339 

Si IHA1207: 

1 3939.5 2.2652 24357 3.6 ± 5.6 ... 13.02 ± 2.41 



a Paramctcrs arc obtained by automatic fitting, assuming C/ — 1.0. 

b Component 5 could not be fit for a continuum that gives physical solution for the other 
components. N V A1243 could be contaminated by a blend or affected by correlated noise in 
this region. 

c It appears that component 2 of C IV A1548 is affected by a poor continuum fit of by 
contamination since C IV A1551 is not detected at the same velocity, so no physical solution 
can be obtained. 



Table 3 

Error Bar on Cf from Continuum Fitting Uncertainty 2 

"(I) (2) (3) (I) (5) (6T~ 

Rb^ 

C f 0.1 0.3 0.5 0.7 0.9 



0.0 



0.1 ... ... ... ... <±0.01 



0.2 










+0.8 
-0.09 


0.3 








< ±0.01 


+0.7 
-0.18 


0.4 








+0.09 

-0.04 


+0.6 
-0.28 


0.5 






< ±0.01 


+0.26 
-0.10 


+ 0.5 

-0.37 


0.6 






+0.04 
-0.03 


+0.4 

-0.15 


+0.4 

-0.47 


0.7 




< ±0.01 


+0.09 
-0.06 


+0.3 
-0.22 


+0.3 

-0.57 


0.8 




+0.03 
-0.02 


+0.17 
-0.10 


+0.2 

-0.29 


+0.2 

-0.67 


0.9 


< ±0.01 


+0.06 
-0.04 


+0.1 
-0.14 


+0.1 
-0.36 


+0.1 

-0.78 


1.0 


+0 

-0.02 


+0 

-0.07 


+0 

-0.18 


+0 

-0.44 


+0 

-0.87 



a Blank entries indicate unphysical combinations of Cf and 
Rb- 

h Rb is the normalized residual intensity of the blue member 
of the doublet, relative to the effective continuum 



